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Abstract 

> ... 

C^i ■ We supplement a very recent paper of R. Crandall concerned with the mul- 

tiprecision computation of several important special functions and numbers. 
I/-) J We show an alternative series representation for the Riemann and Hurwitz 

^. ' zeta functions providing analytic continuation through out the whole complex 

plane. Additionally we demonstrate some series representations for the ini- 
tial Stieltjes constants appearing in the Laurent expansion of the Hurwitz zeta 
function. A particular point of elaboration in these developments is the hy- 
pergeometric form and its equivalents for certain derivatives of the incomplete 
Gamma function. Finally, we evaluate certain integrals including J* RcS=c ^p-ds 
and J* RcS=c -^r-ds, with £ the Riemann zeta function and rj its alternating form. 
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Introduction and statement of results 

Very recently R. Crandall has described series-based algorithms for computing 
multiprecision values of several important functions appearing in analytic number 
theory and the theory of special functions j9]. These functions include the Lerch 
transcendent <J>(z, s, a) = Yl™=o z n /(n+a) s , and many of its special cases, and Epstein 
zeta functions. 

In the following we let B n (x) and E n (x) be the Bernoulli and Euler polynomials, 
respectively, T the Gamma function, and Y(x,y) the incomplete Gamma function 
(e.g., pQ|2j[T3]). An example series representation from [9] (27) is that for the Hurwitz 
zeta function ((s,a), 

as a) - J- V rM(n + Q)] + J- fv ir Bro(a) Am+S-1 (1 1) 

with free parameter A G [0,27r). 

We complement the presentation of [9] with further representation of the Riemann 
zeta function ((s) = £(s, 1) and some results on the Stieltjes constants. The latter 
constants 7 n (a) (e.g., [U [5j El [7]) appear in the Laurent expansion 

((s, a) = + £ K —Tln{a){s - If, (1.2) 

S — 1 *-^ TV. 

n=0 

where 70(a) = —ip(a), the Euler constant 70 = 7 = — ^>(1), and by convention one 
takes 7^ = 7fc(l). Here ip = Y' /Y denotes the digamma function (e.g., [U[2l|T3|). The 
Stieltjes constants may be expressed via the limit formula 



lk{a) = lim 

TV— >oo 



'^ ln fc (n + a) \n k+1 (N + a) 
2-*> n + a k + 1 

n=0 



We let r](s) = (1 — 2 l ~ s )((s) denote the alternating Riemann zeta function, given 
for Re s > by r](s) = Yl™=i ^—^s — • The following series representation provides an 
analytic continuation of the Riemann zeta function to the whole complex plane. 
Proposition 1. Let |A| < it. Then 

^— ' m s 2 *-^ n\ (n + s) 

m=l n=0 v ' 

This representation holds in all of C. In particular, it delivers the special values 
r)(—j) = (—iyEj(0)/2 for j > 0. An alternative expression for the values ^(O) is [T] 
(p. 805) E n (0) = — 2(n + l) _1 (2 n+1 — l)B n+ i for n > 1 in terms of Bernoulli numbers 
B n = B n (0). 

Remark. The values E n (0) = —E n (l) for n > 2 even are zero while Eq(0) = 1. 
The signs of two adjacent values of odd indices are opposite. 

Proposition 1 has many implications. For example, for A = 1 we have the follow- 
ing, with Li s denoting the polylogarithm function. 
Corollary 1. 

n=0 

and 

IcV-Hl + J-l-U.^y-tj^ (1.46) 

We apply the representation (1.1) for a later Proposition. First, we verify some 
properties of the Hurwitz zeta function from (1.1). 
Corollary 2. (a) 

C(s,a)=C(s,a + l) + a~~ s , (1.5) 



(b) 



(c) for integers q > 2, 



and (d) for Re s < 1, 



d a ((s,a) = -s((s + l,a), (1.6) 



Ec(s,-)=0f-l)C(s), (1.7) 

r=l ^ ' 



1 

C(s,a)rfa = 0. (1.8) 

'o 

By induction, (1.6) gives di((s,a) = (—iy(s)j((s + j,a), where (z) n = T(z + 
n)/T(z) is the Pochhammer symbol. 

We let p F q denote the generalized hypergeometric function (e.g., [US [13]). Deriva- 
tives of the incomplete Gamma function with respect to the first argument may be 
expressed in hypergeometric form. An example is given in the following. 
Proposition 2. 

a roo 

T W (s, x) = —T(s, x) = / f- 1 In t e-'dt 
ds J x 

= ^ 2 F 2 (s,s;s + l,s + l;-x) + r(s)[-lnz + ip(s)] + lnxT(s,x). (1.9) 

Later we provide separate discussion for such derivatives. 
Proposition 3. Let Re a > 0. Then (a) 

-7 " tf(a) = 7o(a) - 7 = e-* (-, 1, a) + f) HT^l (1 . 10) 

\e / *— ' ml m 

x ' m=l 

(b) 

T 8 , /|M .C(2) ,* ^ r(0,n + a) ^ (-1)™ S m (a) 

— + 7^a +— 7i a = > ■ > : — , (1.11) 

2 2 *-^ n + a ^-^ ml m 2 

n=0 m=l 



and (c) 



7 3 C(2) /V,C(2)\,, v, ( , C(3), 72(a) 

-6--7^--[y + ^-)^(«) + ryi(«)--3- + ^- 



00 , 1 

^ -3^3(1, 1,1; 2,2,2; -n- a) + -^ 



ra=0 



7 ln(n + a) + y + ^-i + - ln 2 (n + a) 



+£ 



m=l 



-I)™ B m fa) 



m! m 



3 ' 



:i.i2) 



Corollary 3. For A e [0,2vr), 

00 

in r(x) = 7(1 - x) - in \(x - 1) + + Y^ {r[°> K n + x )] - r[o, A(n + 1)]} 



n=0 



+ E 7 i\ f [ g m(a:) - BmjA™- 1 . 
^^ (m — ljml 



m=2 



Alternative forms of the Bernoulli polynomial sums of Proposition 3 and of another 
sum in (1.11) are given in the following. 
Lemma 1. Let Re a > 0. Then (a) 



E ^^ g ^ --7[e- a 2 F 1 (l,a;a + l;e- 1 ) + a(7 + V'(a))]. 



m=l 

In particular, 



mi m 



-y^l = V— = l-ln(e-l). 

m=l n=l 



;i.i3) 



;i.i4) 



(b) 



E 



ra=l 



£ m (a) t m 1 



ml m a 



- [e at 2 Fi(l, a; a + 1; e*) + a( 7 + ln(-t) + ^(a))] . (1.15) 



£ 



ra=l 



Bm(fl) *™ 

m! m 2 



fit 



-/ [e at 2 F 1 (l,a;a + l;e*) + a(7 + ln(-t)+^(a))| — . (1.16) 
a Jo t 



(d) 

^£(^+0) = _l p«- l(1 , a;1 + „;„)*, (1 . 17) 

^^ n + a ojn hiw 

71=0 

Proof of Propositions 

Proposition 1. We apply the integration technique Crandall refers to as "Bernoulli 
splitting" and use the integral representation for Re s > (e.g., p£], p. 807) 

poo fS—1 

T(s) V (s) = / - T —dt. (2.1) 

JO e "+" L 

Splitting the integral at A, we use the generating function 

n=0 

for the integration on [0, A). For the other integration we use a standard integral 
representation for the incomplete Gamma function, together with a geometric series 
expansion, 

poo fs—l °° poo °° / i\m 

/ iri* = E(- 1 ) m / e ' (m+1) ^"^ = E7^r? r[s ' A(m+1)] - (2 - 3) 

JA m=0 JA m=0 ^ ' 

The sum in (1.3) converges uniformly away from poles so provides an analytic con- 
tinuation to the left of Re s > 0. 

For Corollary 1 we recall that for n > [13] (p. 941) 

Tfl 

T(n + l,x)=n\e- x V ^-r, (2.4) 

*— ' m! 

m=0 

and Li 2 (-2) = ^^=i ^ n /™ 2 - 

Remarks. We have supplied details for Crandall's Algorithm 1 for the rj function. 
As far as the values E n (0), their asymptotic form is a case of the following [TT] 
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(24.11.6), 



(-1) 



L(n+l)/2j ^ 



n+1 



£f x ) _+ <( ^W > H eVen > 

4(n!) " J cos(7rx) , n odd, 



and one has the bounds for < x < 1/2 [11J (24.9.5) 

4(2n-l)!2 2n -l 



n 2n 2 2n - 2 

Similar treatment can be made of 

£S-lg— (a— l)t 



> (-If E 2n . 1 (x) > 0. 



e 4 + l 



-dt = 2"T(s) 



c * 



CM, 



+ 1 



r«E 



-l) n 



n=0 



n + a) 



(2.5) 



where Re a > and Re s > 0. This yields, for |A| < n 

a + 



2-T(s 



^ S ^)- ( S ' 2 



/ ^ (m + a) s 2^ n\ (n + s) 

/ J m=l v ' n=0 v ' 

(2.6) 



again providing analytic continuation to all of C. Here the polar part is absent, and 
this representation could be used to develop expressions for the differences of Stieltjes 
constants 7^(0/2) — 7fc[(a + l)/2]. 

Corollary 2. For part (a), use the property B m (a + 1) = B m (a) + ma m ~ l and the 
sum (by pj, p. 941) 

„m \m+s 

_ x a a 



£(-i)' 



a- s [r(s,aA)-T(s)]. 



(2.7) 



171=0 



ml [m + s) 

For part (b), use the derivative dT(a,x)/dx = —x a ~ 1 e~ x and the functional rela- 
tions aT(a, x) = T(a + l,x) — x a e~ x and T(s) = T(s + l)/s. For part (c), use 
YT=\ B m (§) = (<? 1_m - l)B m and (-l) m B m (l) = B m . The following decomposi- 
tion, wherein the generating function (2.23) is employed, can then be used to verify 



the stated property. 



qX(s) 



q s roc y s-l 



T(s) 



e"-l 



dy 



T(s) 



x/q y 8 ' 1 a a r y 3 ' 1 * 

d v+ / -7, — T d v 







ev - 1 



X/q 



e* - 1 



R \m+s— 1 „s 

l-m -°m ^ 9 

r(s) ^— i to! (m + s — 1) r(s) 



■E' 

m=0 



E, 

n=0 " /A /9 



00 „-(n+l)qy ( „qy 



% 



i-i 



e^ - 1 



-dj/. (2.8) 



For part (d), we first note that L B m (a)da = for all to > 1, giving from (1.1) 
r(s) / ({s,a)da = 22 I v ,' v N _ ;J da -f 



n=0 



(ti + a)' 



s-1 



(2.9) 



Integrating by parts, 



da 



1 Z" 1 ^/ w mM 1 

7 — / 1 (s, Am + a)\ —-. r — - 

(s-1) J \da (n + a) 5 " 1 



do 



(n + a) 1 

= — *— \-\ s - l e~ Xn + A s - 1 e- A(n+1) + n 1_ T(s, An) - (n+ l) 1_ T(s, A(n + 1))1 . 
s — 1 

(2.10) 

Then the sum in (2.9) telescopes to — A s_1 /( s — 1) and (1.8) follows. □ 

Remark. More generally than part (c), we have the reciprocity relation ([?], 

Lemma 1) 



Ec(..f-» 



E< 



+ p qb 



p p 



(2.11) 



Here p > 1 and g > 1 are integers, 6 > 0, and min(p/g, q/p) > b. From (2.11) follows 



;[7], Corollary 3) 



r=l 



E«.(7-' 



1-ra p— 1 
^=0 



l + (£-6)- 
p 



(2.12) 



Proposition 2. First integrating by parts we have 

d r°° 

T (1) (s,x) = —T(s,x) = f-Hnte-'dt 

ds J x 

= -/ lnt—T(s,t)dt= -~^dt + lnxT(s,x). (2.13) 

Jx Ot J x t 



We next recall a hypergeometric form of T(x,y), 



x 



T(a,x) = r(a) iFi(a, 1 + a; —x) 



a 



x" v^ a. (—xV 



" PJ (a+j) r- 
with 1F1 the confluent hypergeometric function. This expression may be integrated 

on a finite interval. For an improper integral extending to infinity we need some 

asymptotic information contained in the following. 

Lemma 2. (Asymptotic form of special p F p functions) As x — > oo (a) 

s 2 
2 F 2 {s, s; s + 1, s + 1; -x) ~ e" x — + arVr(s)pnz - if>(s)], (2.15) 

and (b) 

3 F 3 {s, s, s; s+1, s+1, s+1; -x) e- x ^+x- s ^-T(s)[\n 2 x-2\nxip(s)+ip 2 (s)+'4j\s)}. 

(2.16) 
Here ip' is the trigamma function. 

(a) We use a general procedure based upon the Barnes integral representation of 
p F q ([14], Section 2.3). We have 

2 F 2 (s, s;s + l,s+l;-x) = — T(-y)T(y + l)g{y)x y dy, (2.17) 

2m Jin 



where the path of integration is a Barnes contour, indented to the left of the origin 
but staying to the right of — s, and 

*+™- &:ux? -i£v- (2 - i8) 

The contour can be thought of as closed in the right half plane, over a semicircle 
of infinite radius. We then move the contour to the left of y = —s, picking up the 
residue there. We find 

S 2 Res y= _ s T -}^% = ^r(s)[lnx-iP(s)}. (2.19) 

(y + s) z x s 

This gives the algebraic part of the asymptotic form of the 2-F2 function, that is the 
leading portion. The higher order terms come from the exponential expansion of these 
functions, and they are infinite in number. The latter expansion may be determined 
according to ([2], Section 2.3, p. 57). 
(b) proceeds similarly, where now 

s 3 Res„ T }~ V ^l = — r(g)[ln 2 x -2\nxip(s) + tfj 2 (s) + ^'(s)}. (2.20) 

(y + s)* 2x s 

We may integrate in (2.11) from x to b using (2.12), 

b T(s t) x s b s 

V dt = — 2 F 2 (s, s; s+1, s+1; -x) — - 2 F 2 (s, s; s+1, s+1; -b) +T(s) (In b-\nx). 
t s z s z 

(2.21) 

Then with the Lemma in hand, taking b — > 00 completes the Proposition. □ 

Proposition 3. We take A = 1 in (1.1) and, making use of Proposition 2 and 

related considerations, expand the function 

TV W ^ l ^ r(s,n + q) ^ (-ir B m (a) 

T{s)({s,a) 7 = y^~i — i — \^ + l^ \~ 7 1 TT 2 ' 22 

s — 1 *-^ (n + a) s *— ' ml [m + s — 1) 

n=0 v ' m=\ K ' 
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x 



about s — 1, and apply the denning series (1.2). 

Corollary 3. Maintaining the free parameter A, this follows from — L \y+ip(a)]da = 
7(1 — x) — lnT(x), the corresponding expression similar to (1.10), 

-lnA- 7 -^(a) = e-^fl,l,a")+f; ( " 1 | m ^ (a) A- (2.23) 



m=l 



and the property J B m (a)da = B m+ i(a) / (m +1). 

Lemma 1. For (a)-(c), repeatedly integrate the generating function 



.,«.-. 



y^^M z n = ^__ 1 N<27r . (224) 

^ n\ e z - 1 " v y 

n=l 

For (d), write 

^ r(0,n + q) = yj_ r & - t dt_ = y J_ T e -(n +a)v dv 

Z^ n + a Z^ n + a J n+a t f^n + a], 

Then interchange summation and integration and perform another change of variable. 

□ 
Remark. As T(0, x) = — Ei(— x), where Ei is the exponential integral having many 
integral representations, there are a multitude of ways of obtaining (1.17). 

Discussion: Derivatives of the incomplete Gamma function 

Derivatives of the incomplete Gamma function are important for the computa- 
tional methods of [9], and they provide a starting point for many useful integrals. 
Geddes et al. [12] investigated these derivatives, introducing a function T(m,a,z), 
with T(a, z) = zT(2, a, z), m > 1 an integer and initially \z\ < 1, such that 

r^(a, x) = xT(3, a, x) + In 2 T(a,x), (3.1a) 

11 



r^(o, x) = In 2 x T(a,x) + 2x[lnx T(3, a,x) + T(4, a, x)], (3.16) 

and more generally, with P- = il/(i — j)\, 

JmrV \ rn—1 

' ' " ■ = T (m) (a, x) = ln m x T(a, x)+mxJ2 P T~ X ln m " i_1 x T(3 + z, a, ar). (3.2) 



rfa r ' 



i=0 



The function T satisfies the derivative relations 



dT(m, a, z) 

da 



In zT(m, a, z) + (m — l)T(m + 1, a, z), 



and 



dT(m, a, z) 
dz 



■-\T{m — 1, a, z) + T"(m, a, z)]. 



(3.3a) 



(3.36) 



One form of T is [12] (37) 
T(m, a, z) = — Res s= _i 



s + 1 



171—1 



T(a - 1 - sV 



i=0 



■IV 



j„a+i— 1 



«! — a — I) 



~ v (3-4) 



where a is neither zero or a negative integer. This can be obtained by expressing 
T in terms of the Meijer G- function and then using a contour integral for the latter 
function. Using [12] (38) and the relation (o)j/(o + 1), = ajia + i), as we have 
previously done, we may write 

m-2 



\m— 1 



T(m, a, z) 



(m-2)! \dt 



d 



[T(a - ty- 1 ] 



t=o 



v a-l 



+ ("1 



,m— 1 



- m _iF m _i(a, a, . . . , a; a + 1, a + 1, . . . , a + 1; -z). (3.5) 



Proposition 2 and like considerations are in accord with this result. 
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Certain zeta function and other integrals 

Elsewhere we have considered second moment integrals of the Riemann zeta and 
other functions [8J. The following concerns the integrals 



OO 



and 

r / x _ f°° Li c+i t(x) . f Li,(x) , , 

</-oo u ~ "' Jr cS =c a 

Proposition 4. We have 1(c) = —n for < c < 1, and 1(c) — 7T for c > 1, 
7 a (c = 0) = 27r and / a (c) = vr for c > 0, and iz,(c) = 7rx for |x| < 1, x 6 R and c > 
and Il(c) = — nx(l + x)/(l — x) for c < 0. 

Proof. These results may be developed with the aid of Perron's formula [3] (p. 
245), one form of which is [10J (p. 147) 

J_ f (lY d ± = 6 ( x -n), (4.4) 

where 9(x) is the step function taking values 1, 1/2, and as x > 0, x — 1, and 
x < 0, respectively. The line of integration may be moved by including the poles 
of the integrand. This may be seen by integrating over a rectangular contour with 
corners at c ± iT and b ± iT. For the Hurwitz zeta function one has the growth 
estimates [15] (p. 276), with s = a + it and 1/2 > 5 > 0, 

C(s,a) = 0(l) for a > 1 + 5, 
13 



= 0(|t| 1/2 - CT ) for a < 8, 
= 0(|£| 1/2 ) ior5<a<l-5, 
0{\t\ x ~° In |t|) for 1 - 5 < a < 1 + 5, 
= 0(|t| 1/2 ln|t|) for -5<a<S. 



(4.5) 



Further inequalities for ( a (l — s) = S«°=i e 2nma /n 1 ~ s are also given in [15] . Cases of 
these may be used to show that the 'top' and 'bottom' contributions of the rectangular 
contour integration vanish as T — > oo. Then the Cauchy residue theorem is applied. 
Example estimations along horizontal line segments include the following, wherein we 
note the independence of the parameter a, and c\ and c 2 are positive constants. 



C{<r±iT) 



d. 



a 



< 



IC(*±iT)| 



da 



-5 



1 



c\ 



c\ 



<-/ \((a±iT)\da<— I |T| 1/2 ln \T\da = 25— ^r In \T\ 



and 



1+5 

1-5 

1+5 



C(a±iT) 



d 



(7 



< 



1+5 



|C(a±zT)| 



da 



! -„ « 

1+5 



< 



- \((<T±iT)\da<^ \T\ l -°]n\T\da 

L Jl-S J J 1-5 



< 



c 2 lnT 



1+5 



T\ 1 - a da = —{T s -T- s ) 



T 71-5 

The interval [c, b] can then be appropriately decomposed with such contributions. 

The relevant residues are the following. For rj(s)/s, 1/2 at s = 0. For ((s)/s, 
—1/2 at s — 0, being a case of C(0,a) = 1/2 — a, and 1 at s = 1. For Li s (a;)/s, 
x/(l — x) at s = 0. 
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So for instance, for (4.1) with c > 1, the sum on n in (4.4) may be performed. For 

c < 1 the contribution from the pole at s — 1 is subtracted. 

Similarly for (4.3), when c < 0, the stated formula obtains from 2ix (| — jzr)- 
Remark. It appears the integrals J^ It^J dt for p > 1 and < c < 1 take only 

the values —tt or — 2n. 
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